Can the -n + Xd resonance structures be D*D* and D x D molecules? 
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We use QCD sum rules to study the recently observed resonance- like structures in the n + Xci mass 
distribution, £+(4050) and Z+(4250), considered as D*+ D*° and D+D° + D + &1 molecules with the 
quantum number J p — + and J p = 1~ respectively. We consider the contributions of condensates 
up to dimension eight and work at leading order in a s . We obtain mo^D* = (4.15 ± 0.12) GeV, 
around 100 MeV above the D* D* threshold, and m D% p = (4.19 ± 0.22) GeV, around 100 MeV 
below the DiD threshold. We conclude that the D* + D*° state is probably a virtual state that is 
not related with the Z + (4050) resonance-like structure. In the case of the D\D molecular state, 
considering the errors, its mass is consistent with both Z + (4050) and Z + (4250) resonance-like 
structures. Therefore, we conclude that no definite conclusion can be drawn for this state from the 
present analysis. 
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The recent discovery of several missing states and a number of unexpected charmonium like resonances 
in B-factories has revitalized the interest in the espectroscopy of the charmonium states. There is growing 
evidence that at least some of these new states are non conventional cc states, such as mesonic molecules, 
tetraquarks, and/or hybrid mesons. Among these new mesons, some have their masses very close to 
the meson- meson threshold like the A (3872) [l[ and the Z + (4430) Q. Of special importance is the 
appearance of the Z + (4430), observed in the ir + ip' mass spectrum produced in the B° — > K ir + ip' 
decays. Being a charged state it can not be described as ordinary cc meson. Its nature is completely 
open, but an intriguing possibility is the interpretation as tetraquark state or molecular state 0, 0, Q ■ 

The Z + (4430) observation motivated studies of other B° — > K~tt + (cc) decays. In particular, the Belle 
Collaboration has recently reported the observation of two resonance-like structures in the it + Xci mass 
distribution The significance of each of the n + Xci structures exceeds 5a and, if they are interpreted as 
meson states, their minimal quark content must be ccud. They were called Z 1 t ~(4050) and 2^(4250), and 
their masses and widths are M x = (4051 ± 14+J?) MeV, Y x = 82t^±g MeV, M 2 = (4248lgl 1 f^) MeV, 
T 2 = 177+J4+316 MeV. 

There are already theoretical interpretations for these structures as tetraquark states with J p = 1~ 
[3] and as molecular D*D* state with J p = + Q. In this work, due to the closeness of the (4050) 
and Z+(4250) masses to the £)*(2010VD*(2010) and D 1 (2420)5(1865) thresholds respectively, we use 
the QCD sum rules (QCDSR) [1, [l(| EH, to study the two-point functions of the D*D* molecule with 
J p = + , and the D\D molecule with J p = 1~, to see if they can be interpreted as the new observed 
resonances structures ^(4050) and 2^(4250) respectively. Since they were observed in the n + Xci 
channel, the only quantum numbers that are known about them are I G = 1-. 

In previous calculations, the QCDSR approach was used to study the A (3872) considered as a diquark- 
antidiquark state [l2j and as a D*D molecular state (HI, the Z + (4430) meson, considered as a D*D\ 
molecular state @ and as tetraquark states [13] , and the Y mesons considered as molecular and tetraquark 
states [IH . In some cases a very good agreement with the experimental mass was obtained. The QCDSR 
approach was also used to study the existence of a D S D* molecule with J p = 1 + , that would decay into 
J/tpK* — > J/i/jKtt and, therefore, could be easily reconstructed [131 ]. 

Considering the Z+(4050) resonance structure as a D* + D*° molecule with I G J P = 1~0 + , a possible 
current describing such state is given by: 

3 = {daJ l iC a ){cbj' 1 u b ) , (1) 
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where a and b are color indices. 

The sum rule is constructed from the two-point correlation function: 



n(g)=i / d A x e*-*(0\T\j(x)j1(0)]\0). (2) 



On the OPE side, we work at leading order in a s in the operators and consider the contributions 
from condensates up to dimension eight. The correlation function in the OPE side can be written as a 
dispersion relation: 

roo ope, \ 

IT^V) = / ds P - W ( (3) 

J4m% s ~ 1 

where p OPE {s) is given by the imaginary part of the correlation function: np OPE (s) = lm[H OPE (s)]. 
In the phenomenological side, we write a dispersion relation to the correlation function in Eq. 

/ n phen( \ 
ds P - Y + • • • , (4) 
s - q l 

where p phen (s) is the spectral density and the dots represent subtraction terms. The spectral density is 
described, as usual, as a single sharp pole representing the lowest resonance plus a smooth continuum 
representing higher mass states: 

pPhen{s) = x 2 5(s-m 2 D * D ,) + p c ° nt (s), (5) 

where A gives the coupling of the current to the scalar meson D*D*: 

(0\j\D*D*) = A. (6) 

For simplicity, it is assumed that the continuum contribution to the spectral density, p cont (s) in Eq. 
vanishes bellow a certain continuum threshold sq. Above this threshold, it is assumed to be given by the 
result obtained with the OPE. Therefore, one uses the ansatz [16[ 

p cont (s)=p OPE (s)e(s- So ), (7) 

After making a Borel transform to both sides of the sum rule, and transferring the continuum contri- 
bution to the OPE side, the sum rules for the scalar meson Z^, considered as a scalar D*D* molecule, 
up to dimension-eight condensates, using factorization hypothesis, can be written as: 

x 2 e -,nl, D ,/M 2 = r° dg e -s/M* p OPE^ + ^^(M 2 ) , (8) 
Jim 2 c 

where 

p OPE {s) =pP ert ( s ) + p(^( s ) + p ( G2 )(s) + p mlx (s) + p^ 2 (s) , (9) 

with 

P pert W = ^a / J / pl-a-f3)[(a + f 3)ml-af3s}\ 
„« (8)= _3^M "J* | | [(( , + (3)m ,_ Q<i 

p ^ ){s) ^n^p_ J* da j m _ a _ (3)[ia + (3)m 2_ af3sl 

Ctmax 

mix, \ 3m c m^(qq) f da 2 

a mi n 

P { " q) \s) = ^-^1-^1/8, (10) 
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where the integration limits are given by a min = (1 — y/1 — 4m 2 /s)/2, a max = (1 + y/1 — 4m 2 /s)/2, 
fimin — am^/(sa — m^), and we have used (qga.Gq) — m 2 (qq). We have neglected the contribution of 
the dimension-six condensate (g 3 G 3 ), since it is assumed to be suppressed by the loop factor 1/167T 2 . For 
completeness we have also included a part of the dimension-8 condensate contributions 

n ™*< 9 - g)(M2) = _ m 2 c m 2 (qq) 2 j 1 ^ e ^g^ 
^ Jo 

One should note that a complete evaluation of the dimension-8 condensate contributions require more 
involved analysis including a nontrivial choice of the factorization assumption basis [l7j , which is beyond 
the scope of this calculation. 
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FIG. 1: The OPE convergence for the D* D* molecule in the region 2.0 < M 2 < 3.2 GeV 2 for y/sE = 4.6 GeV. 
Comparing to the total contribution, we plot the relative contributions starting with perturbative contribution 
(dashed line), and each other line represents the relative contribution after adding of one extra condensate in 
the expansion: + {qq} (dotted line), + {g 2 G 2 ) (long-dashed line), + m%{qq) (dot-dashed line), + (qq) 2 (line with 
circles), + m 2 (qq) 2 (solid line). 

For a consistent comparison with the results obtained for the other molecular states using the QCDSR 
approach, we have considered here the same values used for the quark masses and condensates as in 
refs. d, [ll GJ, m c (m c ) = (1.23 ± 0.05) GeV, (qq) = -(0.23 ± 0.03) 3 GeV 3 , mg = 0.8 GeV 2 , 

( g 2 G 2 ) = 0.88 GeV 4 , where g = ^4^7. 

To determine the Borel window, we analyse the OPE convergence and the pole contribution: the 
minimum value of the Borel mass is fixed by considering the convergence of the OPE, and the maximum 
value of the Borel mass is determined by imposing that the pole contribution must be bigger than the 
continuum contribution. To fix the continuum threshold range we extract the mass from the sum rule, 
for a given so, and accept such value of so if the obtained mass is in the range 0.4 GeV to 0.6 GeV smaller 
than y/~so~. Using these criteria, we evaluate the sum rules in the Borel range 2.0 < M 2 < 3.5 GeV 2 , and 
in the so range 4.5 < y/so < 4.7 GeV. 

From Fig. [T] we see that for M 2 > 2.5 GeV 2 the contribution of the dimension-8 condensate is less than 
20% of the sum of the other contributions. Using this fact as a criterion to establish a reasonable OPE 
convergence, we fix the lower value of M 2 in the sum rule window as M^ in = 2.5 GeV 2 . 

The comparison between pole and continuum contributions for ^/sq = 4.6 GeV is shown in Fig. [2j 
From this figure we see that the pole contribution is bigger than the continuum for M 2 < 2.9 GeV 2 . 
The maximum value of M 2 for which this constraint is satisfied depends on the value of Sq. The same 
analysis for the other values of the continuum threshold gives M 2 < 2.75 GeV 2 for ^/sq = 4.5 GeV and 
M 2 < 3.1 GeV 2 for ^/sq = 4.7 GeV. In our numerical analysis, we shall then consider the range of M 2 
values from 2.5 GeV 2 until the one allowed by the pole dominance criterion given above. 

To extract the mass ttid+d* we take the derivative of Eq. ([8]) with respect to 1/M 2 , and divide the 
result by Eq. ©. 



a(l - a)M 2 



(11) 
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FIG. 2: The dashed line shows the relative pole contribution (the pole contribution divided by the total, pole 
plus continuum, contribution) and the solid line shows the relative continuum contribution for ^/sq = 4.6 GeV. 
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FIG. 3: The D* D* meson mass as a function of the sum rule parameter (M 2 ) for y/Jo = 4.4 GeV (long-dashed 
line), y^so = 4.5 GeV (dotted line), y/so = 4.6 GeV (solid line) and y/so = 4.7 GeV (dot-dashed line). The crosses 
indicate the upper and lower limits in the Borel region. 

In Fig. [21 we show the D*D* meson mass, for different values of y/so, in the relevant sum rule window, 
with the upper and lower validity limits indicated. From this figure we see that the results are very stable 
as a function of M 2 . 

Using the Borel window, for each value of So, to evaluate the mass of the D*D* meson and then 
varying the value of the continuum threshold in the range 4.5 < y/so < 4.7 GeV, we get mo'D* — 
(4.15 ±0.05) GeV. 

Because of the complex spectrum of the exotic states, some times lower continuum threshold values 
are favorable in order to completely climincte the continuum above the resonance state. Therefore, in 
Fig. [3] we also include the result for the D*D* meson mass for y/s^ = 4.4 GeV. We see that we get a 
very narrow Borel window, and for values of the continuum threshold smaller than 4.4 GeV there is no 
allowed Borel window. Considering then the continuum threshold in the range 4.4 < ^/Iq < 4.7 GeV, we 
get mu-D* = (4.13 ±0.07) GeV. 

To check the dependence of our results with the value of the charm quark mass, we fix ^/sq = 4.6 GeV 
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and vary the charm quark mass in the range m c — (1.23 ± 0.05) GeV. Using 2.5 < M 2 < 2.9 GeV 2 we 
get m D , D , = (4.15 ±0.07) GeV. 

Up to now we have taken the values of the quark-gluon mixed condensate and the gluon condensate 
without allowing any uncertainties. While from Fig. Q] we can see that a change in the gluon condensate 
value has little effect in our results, this is not the case for the quark-gluon mixed condensate. Allowing 
777.Q to vary in the range m 2 , = (0.8 ± 0.1) GeV 2 and fixing ^/sq = 4.6 GeV we get mo'D* = (4.15 ± 
0.06) GeV. Finally assuming a possible violation of the factorization hypothesis, one should multiply 
(qq) 2 in Eqs. (fTU)) and (fTl]) by a factor K. Using K = 2, which means a violation of the factorization 
hypothesis by a factor 2, and y/so = 4.6 GeV we get mwo* = (4.12 ±0.02) GeV. Therefore, taking into 
account the uncertainties in the QCD parameters as discussed above we arrive at 

m D . D . = (4.13 ± 0.07±g : gg t° ° l t° ol) GeV, (12) 

where the first, second, third and forth errors come from the uncertainties in sq, m c , TOq and the 
factorization hypothesis respectively. Adding the errors in quadrature we finally arrive at 

m D , D , = (4.15 ±0.12) GeV, (13) 

where the central value is around 130 MeV above the D*D* (4020) threshold, indicating the existence 
of repulsive interactions between the two D* mesons. Strong interactions effects might lead to repulsive 
interactions that could result in a virtual state above the threshold. Therefore, this structure may or may 
not indicate a resonance. However, considering the errors, it is compatible with the observed (4050) 
resonance mass. 

One can also deduce, from Eq. (J8]) the parameter A defined in Eq. (JBJ). We get: 

A = (4.20l° Jt ±°o-J 6 o ± o1I t°om) x 10- 2 GeV 5 , (14) 

where the first, second, third and forth errors come from the uncertainties in so, m Cl to 2 , and the 
factorization hypothesis respectively. Adding the errors in quadrature we finally arrive at A = (4.20 ± 
0.96) x 10~ 2 GeV 5 . 

Since to obtain the mass we have taken the derivative of the sum rule in Eq. ([8]), it is important also 
to check if the convergence of the OPE and the pole contribution dominance are also satisfied for the 
derivative sum rule. 
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FIG. 4: Same as Fig. 1 for the derivative of Eq. (]8}. 

From Fig. |4] we see that the OPE convergence is even better as from Fig. [T] Therefore, it is correct to 
fix the lower value of M 2 from the convergence of the original sum rule in Eq. ([5]). Regarding the pole 
contribution, we show in Table I the values of M 2 for which the pole contribution is 50% of the total 
contribution, for each value of y/sQ. 
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Table I: Upper limits in the Borel window for the D* D* molecule obtained from the derivative sum rule. 



( GeV) 


ML,. (GeV 2 ) 


4.5 


2.42 


4.6 


2.57 


4.7 


2.72 



From Table I we see that if we impose that both, the original sum rule and the derivative sum rule, 
should satisfy the OPE convergence and the pole dominance criteria, the Borel window exists only for 
■y/so > 4.6 GeV. Therefore, the result for the mass of the D*D* molecule would be even bigger than the 
result in Eq. (|13p . This fact strongly support our interpretation that the D* + D*° state is a virtual state 
and, probably, is not related with the recently observed Z^(4050). This result is in disagreement with 
the findings of ref. @ . 

Considering the ^^(4250) resonance structure as a D\D molecule with I J p = a possible 

current describing such state is given by: 

jil = -J= [(da7M75Ca)(cfc7 5 Ub) + {d a "fzc a ){c b "f ^u b )] . (15) 

In this case, the two-point correlation function is given by: 

O^(g) = ijd±x e*-*(0|T[^(s)jt(0)]|0>. (16) 

Since the current in Eq. is not conserved, we can write the correlation function in Eq. (|16p in terms 
of two independent Lorentz structures: 

U, v {q) = -iMq 2 )^ - Mi) + n (g 2 )M^. (17) 

The two invariant functions, IT and IIo, appearing in Eq. (|17p . have respectively the quantum numbers 
of the spin 1 and mesons. Therefore, we choose to work with the Lorentz structure g^ Vl since it gets 
contributions only from the 1~ state. The sum rule for the meson Z^, considered as a vector D\D 
molecule, in the Lorentz structure can also be given by Eq.® with: 



p pert (s) - 



p<*»>(s) = 

P {G ' 2} ( S )- 
p mlx {s) -- 



p 



2 12 7T 6 

Ct; 

3m c (qq 

2 7 7T 4 

(g 2 G 2 ) 

2 n 7T 6 

3m c ml(qq 
m 2 (qq) 2 



da 
a 



%{l-a-l3) [(a + f3)m 2 c - a/3s] 2 . 



1-a 



dp 

da I t( a + ^) TO c 



a(3s] 



m 2 (l- (a + /3) 2 ) (1- 2a -2/8) 



da 



^(2a + 3/3)[(a 







a /3 s] 



a 



[(a + (3)m 2 c - a/3s] 



2 4 



7T^ 



yl - 4m 2 /i 



n" 



(M 2 ) = - 



rn 2 rn^{qq) 2 
2V 



da 



1 



1M 2 



(18) 



In this case, from Fig. [5] we see that we obtain a reasonable OPE convergence for M 2 > 2.4 GeV 2 . 
Therefore, we fix the lower value of M 2 in the sum rule window as M^ in — 2.4 GeV 2 . The OPE 
convergence obtained from the derivative sum rule is better than the OPE convergence from the original 
sum rule. From the derivative sum rule we get a good OPE convergence for M 2 > 2.2 GeV 2 . 

The upper limits for M 2 for each value of y/so are given in Table II, for the original sum rule and for 
the derivative sum rule. 
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FIG. 5: Same as Fig. 1 for the DiD molecule. 



Table II: Upper limits in the Borel window for D\D molecule with J p — 1 



( GeV) 


Ml. ax (GeV 2 ) (sum rule in Eq.©) 


^maxi GeV 2 ) (derivative sum rule) 


4.5 


2.56 


2.25 


4.6 


2.73 


2.42 


4.7 


2.90 


2.58 



Again we see that the upper limits in the Borel window imposed by the derivative sum rule are smaller 
than the ones obtined with the original sum rule, and this would restrict the range of values allowed for 
the continuum threshold. However, we will allow a small violation in the 50% pole contribution criterion 
in the derivative sum rule, and we will work in the Borel window allowed by the original sum rule. 

In the case of the D\D molecule we get a worse Borel stability than for the D*D*, in the allowed sum 
rule window, as a function of M 2 , as can be seen by Fig. [BJ 
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FIG. 6: The D\D meson mass as a function of the sum rule parameter (M 2 ) for ^/~sq = 4.5 GeV (dotted line), 
y/so = 4.6 GeV (solid line) and y / 5o' = 4.7 GeV (dot-dashed line). The crosses indicate the upper and lower limits 
in the Borel region. 

In Fig. [HI we show the D\D meson mass, for different values of ^/so, in the relevant sum rule window, 
with the upper validity limits indicated. From this figure we see that the results are not very stable as a 
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function of M 2 . 

Using the value of the continuum threshold in the range 4.5 < y/so < 4.7 GeV, and varying m c , ttiq 
and K as discussed above, we get 

™ — (A in+ - 12 + - 05 + - 01 +0-28\ p „ v /-|Q\ 

m DlD - (4.10_ 008 _ 05 _ 09 _ 02 ) GeV, (19) 

where the first, second, third and forth errors come from the uncertainties in sq, m c , TOq and the 
factorization hypothesis respectively. Adding the errors in quadrature we finally arrive at 

m DlD = (4.19 ±0.22) GeV, (20) 

where the central value is around 100 MeV below the Di_D(4285) threshold, and around 60 MeV smaller 
than the mass of the 2^(4250) resonance structure. Therefore, in this case, there is an atractive interac- 
tion between the mesons D\ and D which can lead to the molecular state discussed above. Considering 
the uncertainties in Eq. (gHJ), and the width of the (4250) resonance structure: T 2 = m^lt^lt MeV ' 
it seems to us that it is possible to describe this structure as a D\D molecular state with I G J P = 
quantum numbers. However, considering the uncertainties, the result in Eq. I|20p is also compatible with 
the observed .2^(4050) resonance mass. Therefore, no definite conclusion can be drawn for this state 
from the present analysis. 

For the value of the parameter A defined in Eq. ([6]) we get: 

\z 2 - (l.bU_ 30 _ 30 _ 19 _ 10 j x 10 GeV , (21) 

where the first, second, third and forth errors come from the uncertainties in so, m Cl TOq and the 
factorization hypothesis respectively. Adding the errors in quadrature we finally arrive at Xz 2 = (1.88 ± 
0.77) x 10~ 2 GeV 5 . 

In conclusion, we have presented a QCDSR analysis of the two-point function for possible D*D* and 
D\D molecular states with I G J P = 1~0 + and I G J P = respectively. For the D*D* molecule with 

I G J P = 1~0 + wc got a mass around 130 MeV above the D*D* threshold, and around 100 MeV above the 
observed ^-^(4050) mass. In the case of the DiD molecule with I G J P = 1~1~ wc got a mass around 100 
MeV below the D\D threshold, and around 60 MeV smaller than the observed Z^"(4250) mass. In ref. [2(| 
it was found that the inclusion of the width, in the phenomenological side of the sum rule, increases the 
obtained mass for molecular states. This means that the introduction of the width in our calculation, will 
increase the mass of the D*D* and D\D molecules. As a result, the mass of the D\D molecule will be 
closer to the observed 2' + (4250) mass, and the mass of the D*D* molecule will be far from the Z + (4050) 
mass. Therefore, we conclude that the D* + D*° state is probably a virtual state that is not related with 
the (4050) resonance-like structure recently observed by the Belle Coll. @. Considering the fact that 
the D*D* threshold (4020) is so close to the Z+(4050) mass and that the rj'^&So) mass is predicted 
to be around 4050 MeV [p], it is probable that the ^+(4050) is only a threshold effect [Hj]. Another 
possibility is that the D*D* are in a relative p-wave state. Such configuration will lead to a J p = 1~ 
state as the naive s-wave decay of tt + Xc\ would imply. In the case of the D\D state, although from 
the present analysis its mass is compatible with both, the (4050) and the ^(4250) resonance-like 
structures, in ref. [2(| it was shown that considering a width F = 60 MeV and y/fo = (4.55 ± 0.05) GeV 
the mass obtained for the D\D state is todiD = (4.27 ± 0.03) GeV, in an excelent agreement with the 
Z% (4250) mass. Therefore, their conclusion is that it is possible to describe the Z^ ^(4250) resonance 
structure as a D\D molecular state with I G J P = 1~1~ quantum numbers. 
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